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We investigate the dependence of the stickiness parameters = 1/(12^ ) - where 
the Tij are the conventional Baxter parameters - on the solute diameters cij and o~j 
in multicomponent sticky hard sphere (SHS) models for fluid mixtures of mesoscopic 
neutral particles. A variety of simple but realistic interaction potentials, utilized 
in the literature to model short-ranged attractions present in real solutions of col- 
loids or reverse micelles, is reviewed. We consider: i) van der Waals attractions, 
ii) hard-sphere-depletion forces, iii) polymer-coated colloids, iv) solvation effects (in 
particular hydrophobic bonding and attractions between reverse micelles of water- 
in-oil microemulsions). We map each of these potentials onto an equivalent SHS 
model, by requiring the equality of the second virial coefficients. The main finding 
is that, for most of the potentials considered, the size-dependence of tij(T,o~i,aj) 
can be approximated by essentially the same expression, i.e. a simple polynomial 
in the variable <Ji<jj/afj, with coefficients depending on the temperature T, or - for 
depletion interactions - on the packing fraction 770 of the depletant particles. 
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I. INTRODUCTION 

Theoretical investigation of solutions of mesoscopic particles - with sizes within the range 
10 - 10 4 A - such as colloids, micelles and globular proteins, is more problematic than the 
study of fluids with atomic or simple molecular constituents - with sizes within the range 
1-10 A The main difficulties are due to the large difference between solute 

and solvent molecular sizes, as well to the possible presence of high electric charges and 
large charge-asymmetries. Treating mixtures of macroions and microions, with strong long- 
ranged Coulombic forces, represents a challenge for the most typical methods of the modern 
statistical-mechanical theory of fluids, namely Monte Carlo (MC) or molecular dynamics 
(MD) computer simulations and integral equations (IE) based on the Ornstein-Zernike equa- 



tion coupled with approximate "closures " 



5j. Large size-asymmetries entail very different 



time scales in MD simulations and may lead to ergodicity problems both in MC and MD 
calculations. Moreover, large size-differences imply several difficulties even when using IE 
theories. 

For simplicity, the present paper will be restricted to fluids of neutral particles with spher- 
ically symmetric interactions, neglecting all Coulombic forces due to net electric charges. 
Starting from a fluid mixture with one or more solute species (big particles, or macropar- 
ticles) and one "solvent" species (much smaller molecules or microparticles, which might 
be either a true solvent or polymer coils, smaller colloidal particles, etc.), we will adopt an 
effective fluid approach, which eliminates all large size- asymmetries by averaging out the 
microscopic degrees of freedom corresponding to the solvent . As a consequence, the 
influence of the solvent is incorporated into an effective potential for the interaction between 
big particles, and the initial mixture is reduced to a fluid made up of only solute molecules 
(one or more components). Usually, at the simplest level of description the effective po- 
tential includes, in addition to a steeply repulsive part, a very short-ranged attractive one, 
whose range is a small fraction of the macroparticle size. Recall that a force is said to be 
"short-ranged" if it derives from a potential <pij (r) which vanishes as r~ n with n > 4 when 
r — > oo 0, 9|; the force —d<ftij/dr then decays as r~^ n+1 \ This definition of short-ranged 
potentials is clearly related to the second virial coefficient B 2i ij, which is a central quantity 
in our paper: when the forces are short-ranged in the above-mentioned sense, the integral 
which defines -82,17 ( see ^q- © below) is finite, whereas it diverges for long-ranged interac- 
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tions, i.e. when r < 3. Note that the definition of short-ranged forces is not unique in the 



literature. For instance, in Hirschfelder's classical reference book [10j short-range forces are 
the "valence or chemical forces" , arising from overlap of electron clouds at very short inter- 
molecular separations. The potential of such repulsive, and often highly directional, forces 
varies exponentially with the distance r. On the other hand, all potentials proportional to 
inverse powers of r are called "long-ranged" by Hirschfelder [111 ]. 

Once a reasonable approximation to the effective potential is known, it could be em- 
ployed in both computer simulations or IE calculations. Unfortunately, IEs can be solved 
analytically only in very specific cases, for some potentials and within particular "closures" 
jsj]. The simplest model with both repulsion and attraction which is analytically tractable 
refers to a fluid made up of hard spheres (HS) with an infinitely narrow and infinitely deep 
attractive tail. This hig hly idealized model of adhesive or sticky hard spheres (SHS) was 



proposed by Baxter 



approximation 12 



s mg nl 



13 



and admits an analytical solution within the Percus-Yevick ( 
14j . Notwithstanding its crudeness and known shortcomings 



the SHS model is not a purely academic exercise. In fact, it has seen continuously growing 
interest in the last two decades, because of its ability to describe semiquantitatively many 
properties of real fluids of neutral spherical particles, such as colloidal suspensions, micelles, 
protein solutions, microemulsions, and systems exhibiting phase transitions of several types 



[see for example Refs. 



HQ, 



17 



and references therein) . Accurate simulation data for 
one-component SHS have recently been reported by Miller and Frenkel 19]. 

Because of the simplicity of the SHS model, it has often been suggested to model poten- 
tials comprising a hard core and short-ranged attractive tail by means of sticky potentials. 
To achieve this one needs to define an appropriate equivalence between the actual interac- 
tion and its sticky representation. This mapping of a generic short-ranged potential onto 
a SHS interaction is usually accomplished by requiring the two different models to have 



equal second virial coefficients 



16 



201 ] . Moreover, when applied to mixtures, this approach 
requires a further step, and this is the main poi nt addressed in the present work. 



In a series of earlier papers 



21 



22 



23 



24 



251 ]. we investigated the multi-component SHS 



model, focusing on its possible application to polydisperse colloidal suspensions, namely to 
mixtures where the number p of components is so large that it can effectively be regarded as 
stemming from a continuous distribution. This is, for instance, the case of size polydispersity, 
where - in the discrete notation - a SHS mixture is fully characterized by two sets of 
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parameters, i.e., the HS diameters {<7j} and the "stickiness" coefficients {Uj = l/(12ry)} 
(Tij are Baxter's parameters); the latter depend on temperature T and the strength of the 
interparticle adhesion. Intuitively, one expects to depend on the diameters <7j and Oj of 
the interacting particles % and j, but it is not easy to specify a priori the correct functional 
form, and in our previous papers we attempted some reasonably motivated choices for such 
a dependence. 

The main purpose of the present paper is to investigate the relationship between stickiness 
coefficients and particles sizes, and thus to get new insights into the possible forms of the 
function = (T, <7j, <jj) , starting from a physically sound basis. To achieve this, we will 
present an overview of the most important short-ranged attractive interactions occurring in 
real solutions of colloids or micelles. In doing this, our claim is not to be fully exhaustive, but 
rather to gather sufficient physical information about the mechanisms which cause short- 
ranged attractive interactions in solutions of mesoscopic particles, and the corresponding 
simplest model potentials used for their representation. 

By considering several different systems - dispersion forces, depletion forces, polymer- 
coated colloids, solvation forces (in particular, hydrophobic interactions, and reverse mi- 
celles in water-in-oil microemulsions) - we have surprisingly found strong similarities among 
the simplest models employed to represent this wide variety of physical phenomena. By 
constructing, for each of the relevant potentials, an equivalent SHS representation, we will 
deduce and compare the corresponding expressions for Uj = (T, <jj, <jj). 

The paper is organized as follows. In Section II we will introduce the basic formalism, 
concerning the second virial coefficient, the Baxter SHS model, and the mapping rule for 
getting the equivalent SHS potential from a given short-ranged attraction. Sections III is 
dedicated to the direct van der Waals interaction, while Sections IV, V and VI survey the 
most important short-ranged attractions that are indirect, i.e. mediated by the solvent. The 
hydrophobic effect and interactions between reverse micelles will be considered in Section VI, 
as particular cases of solvation forces. For each model potential, a reasonable approximation 
to the corresponding = Uj (T, 0{, <jj) will be calculated. Finally, a summary, with a brief 
discussion, and our conclusions will be given in Section VII. 
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II. BASIC FORMALISM 



A. Second virial coefficient 



For a multi-component fluid, the second virial coefficient of the osmotic pressure reads 
B 2 = ■ XiXjB 2 ,ij , where Xi is the molar fraction of species i, and the partial second virial 
coefficient for the i — j interaction is given by 

■y p r+oo 

B 2 ,ij = --J fij(r) dr = -27T J Ur)r 2 dr, (1) 

with 

^(r)=exp[-^(r)]-l (2) 

being the Mayer function, /3 = (fc B T) _1 , k B the Boltzmann constant, and T the absolute 
temperature. 

When the actual potential consists of a hard core plus a short-ranged attractive tail, i.e., 
4>ij (r) = (r) + (f)\j l1 (r), one gets 

75 dHS i retail (n\ 

&2,ij = B 2ij + B 2ij [6) 



2,ij ^2,ij 

CO 

tail 



B% = -27T 



r+oo r+oo 

/ tif{r)r 2 dr = -3 / f^(a tJ x)x 2 dx 



(4) 



where <jj is the HS diameter for particles of species i and we set = (<7j + <Jj) /2 as usual, 
introducing also the shorthands B^ = (27r/3) of- and f-j ll (r) = exp [— (r)] — 1. 



Often, the required integration cannot be performed analytically, but if <Plf l {r) is suffi- 



ciently small compared to the thermal energy ksT, then approximate analytical expressions 
may be obtained after expanding the Mayer function flf A (r) in powers of Y = —{3(f)\j l1 (r). 
A numerical estimate of the range of applicability and the maximum relative error A max = 
max |1 — /approx//|, for each of the three simplest approximations, is 



f = e Y -l 



Y < Y < 0.1 A 



max 



/o 



Y + Y 2 /2 0<F<0.6 A max ~5% (5) 
Y + Y 2 /2 + F 3 /6 < Y < 1 A max < 



B. Adhesive hard spheres as a limiting case of square-well model 



Probably, the simplest two-parameter representation of a spherically symmetric inter- 
action with steeply repulsive core and short-ranged attractive tail is the square-well (SW) 
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potential 

+00 < r < (Tij , 



r 



o"« < r < a {j + w i7 - , (6) 



r > cry + Wij , 



with £jj > and w^j being the depth and width of the well, respectively. The corresponding 
partial second virial coefficient reads 

BlJ 3 = {1 - (e*« - 1) [(1 + A.-) 3 - 1] } (7) 



R HS 



1 " 3 (e*« - 1) (a,, + A» + ~Aj.) 



with Ay = Wij/dij > 0. Eq. (JTj) shows that, if the well is narrow (Ay <C 1), -Bf^ can be 
significantly different from -B^- only when the attraction is strong enough (e^ 3> l) . 

Unfortunately, despite the simplicity of the SW model, no satisfactory analytical solution 
of the resulting IEs has been found so far. However, such a solution can be found within 
the Percus-Yevick (PY) approximation for a special limiting case, when the well width Ay 
goes to zero but the depth goes to infinity in such a way that the contribution of the 
attraction to the second virial coefficient remains finite and different from zero (Baxter's 

n 

sticky limit) [12] . The short-ranged attraction becomes a surface adhesion, and the particles 
of the resulting model are thus named adhesive or sticky hard spheres. From Eq. (J7|) one sees 
that Baxter's condition on Bf^ requires the product (e^ e « — l) A^- = to be independent 
of Ay for small A^-, and this leads to the following condition for the SW depth 



^Baxter SW = A . flT ^ ( 1 + ^_ ) (g) 



As previously mentioned, our is simply related to Baxter's original parameter by 

Here, measures the strength of surface adhesiveness or "stickiness" between particles of 
species i and j, and must be an unspecified decreasing function of T. In fact, as T — > 00 one 
must also have —>■ 00, in order to recover the correct HS limit. The SHS models must 
therefore satisfy the high-temperature condition 



lim t tj = 0. (10) 

T — >oo 
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A consequence of Eq. (JSJ) is a very simple expression for the SW Mayer function 

-1 < r < a 



Baxter SW/ 
ij \ 



Uj Vij/Wij CTjj <r< (Tjj + Wij , (11) 

r > Oij + , 



Baxter focused on fy, since this quantity directly determines i?2,tj and, furthermore, the 
coefficients in the cluster expansion of thermodynamic properties and correlation functions 
can be expressed in terms of mult i- dimensional integrals of products of Mayer functions . 



The simple functional form of f i3 ter bW (r) then allows one to calculate analytically many 
quantities of interest. In the "sticky limit" {wij} — > {0}, the Mayer function becomes 

/ y SHS(r) = [9(r- cry) - 1] + % a tJ 5 + (r - o %j ) (12) 

with 9(x) being the Heaviside function (= when x < 0, and = 1 when x > 0) and 5+(x) 



an asymmetric Dirac distribution 26], while the SHS second virial coefficient is simply 



B 2*ij = B 2% (1 - 3tjj) • (13) 



C. Mapping onto equivalent SHS model 

On comparing Eqs. (|T3l) and (jl]), one has 

rjtail 

.eq(tail) _ ■ D 2,ij /., ,\ 

13 " 3S 2 H|' U4j 

and hence the following mapping rule: the parameters Uj of the equivalent SHS model must 
be given by 

i />+oo r+oo 

^ (tail) = 4 / ftfirVdr = / ffifaxWdx. (15) 

This is the main relation used in the remaining part of the paper. The superscript in 
means: this tij yields the SHS potential equivalent to <p\j l1 . 



III. VAN DER WAALS ATTRACTION 



The main direct attraction between two neutral molecules % and j is the van der Waals 
(vdW) interaction, represented by the potential 0^/ W ( r ) = — C^ w r -6 , which is - in general 
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- the sum of three different contributions. For most simple molecules - except the small 
highly polar ones - the vdW attraction is almost exclusively determined by the dispersion 
forces; the latter are in fact the only contribution to the vdW forces if both molecules are 
nonpolar. 



A. Dispersion forces 

The dispersion or London forces are induced- dipole/induced-dipole interactions, whose 
potential is given by the London formula flo| 

C 3 II 

4f v ( r ) = -7? ' c u = 2T+T- a ' ia ' j ' for large r ' ^ 16 ^ 

where f and a[ are, respectively, the ionization energy and polarizability volume for 
molecules of species i. As the name suggests, a[ has the dimensions of volume. It can 
also be written as oiij (47T£ ), where e§ is the permittivity of the vacuum and «j is the polar- 
izability of species i, which increases with increasing molecular size and number of electrons. 
Hence the polarizability volume is proportional to the molecular volume, i.e., a[ oc of. 

This polarizability effect alone can produce considerable molecular attraction, and is 
responsible for the formation of liquid phases from gases of nonpolar substances (argon, hy- 
drogen, nitrogen, etc.). The name "dispersion forces" stems from the fact that the electronic 
oscillations producing the London attraction are also responsible for the dispersion of light. 



B. Hamaker's macroscopic approximation 

Colloids, micelles and globular proteins are mesoscopic particles formed by a very large 
number of polarizable molecules (typically 10 10 in micrometer-sized particles) [l|. As a 
consequence, the total attraction energy between such macroparticles can be obtained by 
pairwise summation of London energies between all molecules of the two interacting bodies. 



Hamaker 



2 ] 



for the energy of interaction of 



271 ] performed an approximate calculation 
two fully macroscopic bodies i and j in a vacuum, with densities pi and pj and occupying 
volumes Vi and Vj. Replacing the discrete distribution of molecules inside each body with a 



continuous one, Hamaker obtained for two spheres of arbitrary size 



0g (r) 



A| \ a i<Tj t (JiOj f n ^ fr 2 -o 2 - 
12 



+ 'to +21n 



r 2 - cr 2 r 2 - L\ \r 2 - L 2 - 



12 



(17) 



where = fa — aj\ /2, and cr^- < r < +oo. Here, = n 2 piPjCij 4J is referred to as 
Hamaker's constant, and has dimensions of energy. As oc a^a'- oc cr?cr|, and pjPj oc 
o"^ 3 cr^" 3 , v4^- is nearly independent of i and j. In the case where all mesoscopic particles are 
made up of the same material but have different diameters (discrete size polydispersity) Afj 
reduces to An = ii 2 p 2 C, which is a property of the material itself. 

Hamaker's macroscopic result has also been applied to mesoscopic particles, with the 
justification that the potential (TlTj) has a scaling property: if r, <7j,<7j are all multiplied by 
a factor 7, the attraction energy remains unaltered, i.e., (77-, 7<7j, 70,-) = (r, <7j, cry). 
Note, however, that Hamaker's formula refers to two spheres in free space, i.e., it neglects 
the screening of London forces due to the suspending medium. 

In the limit r — > +00, a series expansion of Eq. (117j) yields 

^3 

36 ~ 



<t>l ( r ) ~ for r » °<i > > ( 18 ) 



which means that at large distances the two spheres behave, to leading order, like point- 
particles even though the factors af and cr| stem from HS volumes. 

On the other hand, the Hamaker potential is singular at contact, i.e., when r — > cr^-. 
This is due to the approximation of regarding the two spheres as continuous distributions 
of point-particles, neglecting all intermolecular repulsions. The leading divergence is 

4>Vj(r ~-T7f- — 1 for 0<r-a jj «min (7 il (7 i . 19 

J 12 /- crfj 24 a, , r - cry 

This divergence simply means that the continuum picture must break down and molecular 
granularity - with excluded-volume effects - cannot be neglected once the closest distance 
r — cry between the two spherical surfaces becomes very small. 

Such a deep attractive potential would lead to irreversible association or "flocculation" 
of the suspended particle. This effect can be avoided in one of two different ways, namely 
by charge stabilization or steric stabilization. In the first case, some surface chemical groups 
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of the particles become partially ionized in water, and the resulting electrostatic repulsion 
makes close contact impossible. In the second case, stabilization is achieved by grafting 
polymer chains ("hair") to the particle surfaces. Both stabilization mechanisms - exten- 
sively used for colloidal suspensions - imply that the closest approach distance between % 
and j becomes larger than cr^, i.e., of? = cr^- + 5, with 6 > being an additional character- 
istic length. The Hamaker singularity at contact is thus avoided, and the vdW attraction 
may then be treated as a small perturbation, if the effective HS diameter is sufficiently large 
compared to the bare one (in sterically stabilized colloidal suspensions, o~ eff exceeds a typi- 
cally by 10%). Moreover, it is possible to strongly reduce the value of the Hamaker constant 
by "refractive index matching" 

A numerical estimate of the strength of the Hamaker attraction is given - in the one- 
component case, for simplicity - by the quantity 

F max = -/30 H (a + 5) = 3 H (1 + A) ^, (20) 

where A ee S/a, 

H(u) = — + \ + 2 In ( 1 - -V), (21) 



u 2 — 1 U 2 \ U 2 j 

and, from Eq. (Tl8|) . we have defined a Hamaker temperature as 

Th = , (22) 

36^ ' v ; 

which depends on the material which constitutes the particles. In most cases, An lies 
between 10~ 20 and 10~ 19 J, i.e., 2 k B T < A u < 20 k B T, where T = 298.15 K. A typical 

value An = 0.5 x 10~ 20 J (= 10 k B T) yields T H = 100 K, and thus - at room temperature - 



-2 k B T if A = 0.1 , 
-0.6 k B T if A = 0.2 , 
-0.2 k B T if A = 0.3 . 



Using Eq. (f2"U|) for Y m3 ^ together with the criteria in Eq. (j^D, one finds the approximate 
lower bound T mm /K for the applicability, respectively, of the linear, quadratic and cubic 
approximations to the Mayer function, as reported in Table [H 

Thus, whereas the linear approximation works only at high temperatures, the quadratic 
one is already sufficient even at room temperature if A > 0.2. 
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Unfortunately, analytical integration of the expression (ITT]) is not possible, and conse- 
quently no result for can be obtained directly from (f)fj (r). Nevertheless, in order to 
get a rough approximation to tlj , we propose an analytically integrable interpolation of 
the correct behavior of cffj ( r ) at short and large distances, i.e. 



/ \ 3 3 r t 

r — o~. 



1 — exp 



L 



(23) 

where + S < r < +oo, and L acts as a screening length. When 5 ~ 0.1<r, L = 0.108a ~ 5 
yields a satisfactory contact value, i.e., 0^~ mterp (a^ + 5) ~ (a^- + 5). Using the linear 
approximation - valid at high temperatures - one gets 



t 



1 A H 

eq(H— intcrp) ij 0~iO~j 



^ ( ^ + 2e- 5 ' L — + (l + 41 ( A X 2 



L / aj j y L y y a %3 



where Ei(z) = ^-du is the exponential integral. However, since the factors L/afj ~ 
S/a™j refer to big particles, the leading term - at least within the linear approximation - is 

eq(H-interp) _ 3 / 5 \ ( T H a^A 

~ 2 E ' [l) [t -3t) ■ (24) 

C. Polarizable hard spheres. Sutherland model 

Focusing only on the r~ 6 part of the Hamaker potential, which represents the long- 
distance polarizability, one could define a simpler model, corresponding to a mixture of 
mesoscopic HS with dispersion attractions, called polarizable hard spheres (PHS), i.e., 

I'us,., ) +°°> 0<r<<7y 
-Aij afa^/r 6 r > <r t . 



'131 



where the choice 



An = (26) 



ensures the mesoscopic size of the particles. If all particles are made up of the same material 
substance, then 0^ HS (r) = —A afaj/r 6 for r > a^. 

The potential ( 1251) may be regarded as a special case of the Sutherland model, which 
represents rigid spheres which attract one another according to an inverse-power law, i.e., 



r) = — e (a/r) for r > a (e > 0) [lOj. Indeed, one could rewrite it as (r) 



— e^ HS (cTij/r) , with e™ s = A (aiO-j/afA for particles with the same material composition. 
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The strength of this interaction - in the one- component case - is then given by 

T H 100 K 



Knav / 



iPHS 



(J 



T 



T 



after taking Th ~ 100 K. From the results outlined in Eq. dSJ), a linearization of the 
Mayer function makes sense for T > 1000 K. A quadratic approximation is feasible when 
T > 167 K. Finally, the cubic approximation holds for T > 100 K. 

Therefore, in the multi-component case, one can safely adopt the cubic approximation to 



/P HS (r) and perform the integration in Eq. ( 1T51) . obtaining 



, C q(PHS) _ 1 9 T H (uiUj \ 3 



07 



2 



(27) 



In the one-component case this expression reduces to 



f q(PHS) = 12 ^ + 15 



7h 



+ 8 



T 



For Th ~ 100 K, this results - at T « 300 K - in a value t ec *( PHS ) ~ 5.96, which corresponds 
in Baxter's parametrization to 



t = — ~ 0.014, 
12t 

and lies well below the critical temperature of the SHS fluid, r c = 0.1133 ± 0.0005 



IV. EXCLUDED VOLUME DEPLETION FORCES 

In general, the indirect, solvent-mediated, solute-solute interactions depend on both the 



so 



28 



ute-solvent and solvent-solvent ones, and thus may be very difficult to evaluate 



7 



2d 



291 ] . We will now report several very simplified cases. 



Asakura and Oosawa (AO) |30j, and independently Vrij [3l|, first showed that two big 
(colloidal, or solute) particles, immersed in a sea of small particles, feel a mutual attraction 
when their surfaces get closer than the size of the smaller particles (depletion attraction). 
This effect is an indirect attraction originating from the interactions of the two big particles 
with the small ones of the environment, even if these latter consist of, say, hard spheres. In 
mixtures with neutral components, the small particles - hereafter referred to also as depletant 
particles - may correspond, for example, to solvent molecules, non-adsorbing polymer coils, 
or smaller colloidal particles. 
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Upon adding, for instance, polymers to a stable colloidal suspension, the colloidal particles 
tend to aggregate. The polymer-induced depletion forces between the colloidal particles can 
cause formation of colloidal crystals or flocculation. 

In the AO model, originally designed to describe colloid-polymer mixtures, the big-big 
(colloid-colloid) interactions as well as the big-small ones are modeled as excluded-volume 
HS interactions, while the small-small interactions are assumed to be zero (ideal gas approx- 
imation, corresponding to mutually interpenetrable, non- interacting depletant molecules). 
In particular, polymer coils are assumed to have an effective HS diameter equal to twice 
their radius of gyration. 

Consider two big HS of species i and j at distance r, with radii Ri and Rj, in a dilute 
suspension of depletant spheres of species 0, with radius i?o- The solute molecule i produces 
a spherical excluded-volume region of radius = Ri + Rq around itself where the centers 
of the depletant particles cannot penetrate; this is also called the depletion zone. When 
the shortest distance r — (Ri + Rj) = r — Oij between the surfaces of % and j becomes less 
than the diameter cq = 2Rq, the two depletion spheres surrounding i and j overlap and 
the small particles are expelled from the region between the big molecules. This implies 
that the thermal impact forces on the pair i and j from the "outside" are only partially 
compensated by those from the "inside" (see Fig. 8 of Ref. |3l|). The depletion effect is due 
to this unbalanced pressure difference, which pushes the big particles together, resulting in 
a net attraction. From another point of view, the overlapping of excluded volumes implies 
that the total free volume accessible to small particles increases, leading to a gain in the 
system entropy with a consequent decrease of the Gibbs free energy. This trend to decrease 
ree energy produces an effective indirect attraction between the big spheres. AO and Vrij 



30. 



3 if ] evaluated the HS- depletion (HS-depl) potential as 

f+oo < r < (Ty 

-p k B T ^° vcrlap (r) a l3 <r<D l3 (28) 



ij V 1 (J — ' — ^13 

r> D i:i , 



^overlap^ = _ r) 2 [3^.. _ ^ ^ _ f) + (% _ ^ (2g) 



where po is the number density of the depletant molecules, = cr^ + o- , and l^° verlap (r) 
denotes the lens-shaped overlap volume of two spheres with radii a i0 = (a - , + a ) /2 and 
°~oj — (ffj + °"o) /2, at distance r (see Appendix). The attraction increases linearly with 
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temperature and with the concentration of depletant particles. Since the AO model includes 
only HS interactions, the corresponding depletion forces have a purely entropic origin. Fi- 
nally, it should be emphasized that the AO approximation is valid only for dilute suspensions 
of depletant molecules, i.e., at low po values; formally, this last restriction can be removed 
by replacing po by the density of polymer in a large reservoir connected to the system. 

The tail of 0^ S ~ dcpl (r) has a finite range, equal to the diameter a of the depletant 
molecules. The attraction strength can be estimated from 



Y: 



i] ,max 



^HS-depl 



a 



1.1 1 



Vol I 



3 (Ji<7j 1 

2 Oij cr 



(30) 



where r/ = (tt/6) po^o is the packing fraction of the depletant particles. Note that Yij jTaax 
does not depend on temperature, but the attraction strength may be tuned by varying i]q. 

For one-component solutes, Yij >max reduces to Y max = r/ (1 + 1.5/A) with A = cr /a. In 
this case - following again the criteria given in Eq. ([5]) - the upper boundary r^™ ax for the 
applicability of the linear, quadratic and cubic approximations to the Mayer function given 
as a function of A, respectively, is reported in Table [III 

The quadratic approximation result for the equivalent SHS model is 
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where Ay = Co/cij- As remarked, this expression does not depend on T, since the solute- 
solvent interactions are of purely HS character. If Ajj <C 1, then the leading terms are 
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Generalizing from the form of the quadratic approximation for general A^, one expects 
the cubic approximation to yield a result of the form 
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Several other studies of depletion forces, which go beyond the entropic HS approach 
by taking into account more refined representions of t 
interactions, are also available in the recent literature 



ie solute-solvent and solvent-solvent 
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V. POLYMER-COATED COLLOIDS OR HAIRY SPHERES 

If the intermolecular attractive forces are strong enough, a colloidal suspension phase- 
separates, or even flocculates or gels. As explained above, stability against flocculation may 
be ensured by steric or charge stabilization. In steric stabilization, the colloidal molecules 
are coated with grafted polymers - the "hair" - which can prevent particles from coming 
sufficiently close to experience a strong vdW-attraction. 

However, changing the solvent or the temperature may turn the effective interaction 



from repulsion (HS-behavior) to attraction 32 



33 



34j. When sterically-stabilized colloidal 



particles are immersed in a good solvent for the polymer brushes, the solutes behave like HS, 
independently of temperature; this is the case, for example, of silica particles covered with 
a layer of octadecyl chains, when immersed in cycloexane. On the other hand, for each poor 



solvent there exists a Flory's theta-temperature Tg [35j, which is characteristic of the given 
solvent-polymer pair and has the following property: the solute particles behave like HS at 
T > Tg, whereas an attraction occurs at T < Tg. This occurs with e.g. silica particles with 
octadecyl chains, when dispersed in benzene. The term "^-solvent" indicates a poor solvent 
at T = Tg. 

These effects originate from a competition between polymer-solvent and polymer-polymer 
interactions. First, the nature of the solvent influences the polymer conformation. In fact, 
in a good solvent the interactions between polymer elements - monomer units - and adja- 
cent solvent molecules are strongly attractive and thus predominate over possible intra-chain 
polymer attractions. Consequently, the polymer will assume an "extended-chain" configura- 
tion, so as to reduce the number of intra-chain contacts between monomer units. Polymer- 
coated colloidal particles will have fully-extended hair and thus the largest HS diameter 
possible, corresponding to the strongest solute-solute repulsion. 

In a poor solvent, on the other hand, the polymer-solvent attractions are weak. Now it 
is the temperature that determines the solute-solute interaction. At T > Tg the hair will be 
fully-extended, as in good solvents (HS behavior). At low temperatures T < Tg the polymer 
segments find their own environment more satisfying than that provided by the solvent. 
This preference may produce more compact "globular" configurations, in which intra-chain 
polymer-polymer contacts occur more frequently ("curly hair"). In an alternative view, when 
two solute particles are in close contact, a high number of polymer-polymer attractions is 
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favored by the interpenetration of the two polymeric lay ers. 



In the literature on sterically-stabilized colloids [32|, [33|, [34j, the attractive part of the 
potential for one-component hairy hard spheres (HHS) - due to the polymer-polymer inter- 
actions between surface layers of different particles - was described by a SW, with a depth 
proportional to the (maximum) overlap volume of the layers and temperature-dependent in 



361 ] . The SW width equals 



analogy with the Flory-Krigbaum model for polymer segments 
the length of interpenetration of the stabilizing chains, whose maximum value coincides with 
the thickness i of polymeric layer. 

For mixtures, we could consider the most direct extension of the one-component SW 
model. In such a case, using Eq. ( 1T51) . a SW potential could immediately be mapped onto 
a SHS one: 
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When the SW is very narrow (A ? ^ w <C l), one can approximate 

^ (SW) - [exp (/?4 W ) - 1] Af. 
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However, instead of a discontinuous SW model, we prefer to propose a potential with a 
similar but continuous attractive tail of finite range, i.e., 
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where 



F(T) 



(35) 



T >T e 

Te/T-l T <T e 

Here, we call F(T) Flory's temperature-function, and 0™ s (r) is assumed to be proportional 
to the overlap volume between polymeric layers of the two HHS at separation r, with pi 
being a number density proportional to the polymer density in the stabilizing layer. 
Within the linear approximation, one finds for the equivalent SHS model 
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where r]i = \pi (2£) 3 . Since the thickness I is much smaller than the particle sizes, one may 
expect - by analogy with the HS-depletion model - that the cubic approximation reads 
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Note that, since F(T) = when T > Tg, then limr^ +00 F(T) = 0. Thus the form of 
F(T) ensures that the high-temperature condition ( fTUI) is satisfied. 

VI. SOLVATION FORCES. GURNEY-FRIEDMAN MODEL 



An indirect interaction between solute particles may also arise from solvation. To picture 
solvation effects, Gurney 37J] and Frank and Evans [38( introduced the physically intuitive 



concept of cosphere or solvation layer. One supposes that any isolated solute particle is 
surrounded by some region in which the solvent has different properties than the bulk solvent, 
since its structure is markedly affected by the presence of the solute: some of the solvent- 
solvent bonds have been broken by the introduction of the "foreign" particle. Clearly, such 
a region has no well-defined boundary, but Gurney's model assumes that significant effects 
come from only the few solvent molecules that are directly next to the solute particles, i.e., 
in a spherical shell whose thickness 5 is taken to be a few solvent diameters or even the size 
of only one solvent molecule (for water, a molecular diameter of 2.76 A is acceptable). This 
picture was first applied to electrolyte solutions by Friedman and coworkers [39|. In the 
ionic case, however, the previous definition of cosphere, with the same thickness for every 
ionic species, may be too restrictive, since the solvation region may extend even outside the 
cosphere, as occurs for very small ions (Li + , and polyvalent ions such as Mg 2+ , Ca 2+ , etc.). 

When two solute particles i and j approach sufficiently closely for their solvation layers 
to overlap, some of the cosphere solvent is displaced. Furthermore, the overlapping region 
contains solvent molecules which are now affected by the combined force field of two solutes, 
so that its structure might even differ from that of each isolated cosphere. The whole process, 
in which the sum of the cosphere volumes is reduced by overlap and the solvent relaxes to 
its normal bulk state, will be accompanied by a Gibbs free energy change. If the solvent 
molecules in the isolated solvation layers are in a state of lower free energy than those in 
the bulk, the overlap of two cospheres with the consequent expulsion of solvent gives rise to 
a free energy increase, and the resulting contribution to the interaction between two solutes 
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is repulsive. When the solvent molecules in the solvation layers are in a state of higher free 
energy than those in the bulk, the expulsion of solvent from the overlapping region leads 
to a free energy decrease. In this case, both the free energy and the disruption of solvent 
structure are minimized when two solute particles % and j are brought close together, causing 
a net i — j attraction. 

Because of the lack of knowledge about the properties of the solvent in the solvation 
region, it is difficult to construct a detailed and physically sound microscopic model of 
the effects described above. Adopting a heuristic approach, Friedman and coworkers [39] 
proposed that the free energy change accompanying the cosphere overlapping of two HS 
solute particles i and j gives rise to the Gurney potential, defined by 



^ Umey (r) = A,(T,p) 



^overlap ( R . + ^ R . + ^ ^ 



V 



(37) 



Here the Gurney parameter A^ is in general a function of temperature T and pressure p 
and represents the molar free energy of transfer of solvent from the overlapping region of 
the i — j cospheres to the bulk. As previously discussed, < corresponds to attraction. 
Furthermore, vq is the mean volume of a solvent molecule, while the volume of solvent 
returning to the bulk is given by the intersection volume of the cospheres surrounding the 
two solute HS at distance r, namely V overlap (Ri + 5, Rj + 5, r). The free-energy parameters 
Aij were determined numerically by fitting the model to experimental data. 

The close resemblance of </>^ urney (^) to both </>™ s (r) and 0^ S_depl (r) is apparent. By 
analogy one obtains immediately 
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This expression may be applied, in particular, to both cases of solvation interactions - 
hydrophobic bonding and interactions between reverse micelles - whose physical origin will 
be illustrated in the following. 



A. Hydrophobic interaction 

The hydrophobic interaction (or hydrophobic bonding) consists in the tendency shown 
by nonpolar portions of molecules or ions with long nonpolar chains or aromatic groups - 
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for example, surfactants, phospholipids, glycerides, and dyestuffs - to aggregate in aqueous 



so 



40 



utions, partially or completely removing such parts from contact with the solvent 



4=11 ] . This hydrophobic attraction between nonpolar entities, which occurs exclusively in 



water, has - to a large extent - an entropic origin, related to the strong tendency of water 
molecules to form hydrogen bonds and associate |4j . 

The physical mechanisms underlying the solvation forces are rather well understood when 
the solvent is water. Polar molecules or polar groups of a solute feel strong attraction towards 
water molecules, and thus are said to be hydrophilic ("water- loving"). On the other hand, 
nonpolar molecules or nonpolar groups "dislike" water, and are called hydrophobic ("water- 
hating", or "water-avoiding"). 

The hydrophobic effect means that nonpolar particles have an extremely weak solubility 
in water, since inserting one of them - a noble gas atom, a chlorine or oxygen molecule, a 
hydrocarbon molecule, etc. - into water may actually lead to an increase of Gibbs free energy, 
i.e., AG so i u tion > 0. Indeed, the formation of a new cavity requires the breaking of many 
water-water hydrogen bonds with a considerable AG > 0, which cannot be compensated 
by the small AG < provided by the new solute-water vdW interactions (the nonpolar 
solute cannot participate in the formation of hydrogen bonds). Then, in order to get a 
further decrease of G, the water molecules close to the solute reorient themselves, so as to 
create as many hydrogen bonds with adjacent water molecules as possible. The result is the 
formation of a "cage" - or hydration layer - around the solute, with more rigid water-water 
bonds than in the bulk. Such an additional ordering in the solvent, brought about by the 
introduction of a solute molecule, implies a significant entropy decrease, AS < 0, and thus 
a strong positive contribution — TAS to the total AG = AH — TAS of solution at constant 
T and p. This explains why nonpolar particles are hydrophobic: one gets AG so i u tion > 0, 
when the entropic contribution dominates over the entalpic one AH, which is usually small 
and may be positive or negative. On the other hand, at higher temperatures the solubility 
increases, and one may find AG so i u tion < 0. In fact, AS becomes much smaller, because the 
molecular thermal motion struggles more efficiently against the structure formation around 
a nonpolar solute. 

The same hydrophobic effect is responsible for the above-mentioned hydrophobic bonding, 
where nonpolar parts of molecules or ions tend to aggregate. In fact, the solvent molecules 
prefer mutual contacts over those with the "foreign" substance (solute), while the aggre- 
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gation of solutes reduces the total volume of their "cages" , minimizing the loss of entropy. 
The hydrophobic interaction arises when overlap of hydration-layers occurs, and becomes 
increasingly attractive as the distance between two solute particles decreases. 

This phenomenon in acqueous solutions of alcohols was studied by Friedman and Kr- 
ishnan 42J, who used a model potential containing a repulsive term of the r~ 9 form, plus 
an attractive Gurney term given by Eq. (157)1 . representing the overlap between cospheres. 
For the sustances they considered, these authors found values of the Gurney coefficient 
A xx in the range —190 to —60 cal mol -1 . The cosphere thickness 5 was taken to be 2.76 A, 
corresponding to one molecular layer of water. 

Clark et al. 43j investigated the same physical systems with a more refined model, 
including a core potential based upon Lennard- Jones potentials for individual atom-atom 
interactions and again a Gurney term for the hydrophobic attraction. They found a mini- 
mum ~ —0.5 ksT in their potentials of average force, that implies an overall tendency for 
those alcohols to associate when in an aqueous environment. 

Hydrophobic bonding is very important in interface and colloid science. It is often the 
driving force behind the way in which biomolecules organize themselves and it is responsible 
for the formation of micelles, bimolecular layers, and lamellar structures. 



B. Reverse micelles in water-in-oil microemulsions 



Molecules having both hydrophilic and hydrophobic parts are said to be amphiphilic 
( "dual- loving" , in the sense of being both "water-loving" and "water-hating"; from the 
Greek ampt = "on both sides"). An important example is provided by relatively short chain 
molecules with an ionizable or polar (thus hydrophilic) head-group and a nonpolar (thus 
hydrophobic) tail, consisting of one or several flexible hydrocarbon chain(s). Since these 
molecules can significantly lower the surface tension of a solution, they are generically called 
surfactants or surface- active agents US , l^lj]. When immersed in water, the head-group 
may become negatively or positively charged, or it remains polar with no net electric charge. 
Accordingly, the surfactants are classified as anionic, cationic or non-ionic. 

Clearly, it is the hydrophilic head-group that keeps a surfactant solute dissolved in the wa- 
ter. The hydrophobic tail tends to avoid contact with water and to seek, as far as possible, a 
nonaqueous environment. The longer the hydrophobic tail, the poorer the solubility in water 
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and, hence, the greater the tendency of the surfactant to escape from the acqueous solution. 
Consequently, as the solute concentration increases, a phase-separation may occur. Alter- 
natively, the surfactant molecules accumulate at interfaces between water and other liquids 
or gases, or spontaneously self-assemble into supramolecular aggregates of mesoscopic size 
that minimize the number of contacts between water and hydrocarbon tails and maximize 
the number of tail-tail interactions. 

If a surfactant is added to pure water under atmospheric pressure, its molecules first form 
a monolayer film at the water-air interface, with polar heads pointing towards the water and 
tails exposed to the air. Above a certain concentration of surfactant (critical micellization 
concentration), one finds an abrupt change in the properties of the solution; in fact, now 
the solute particles in the bulk begin to form supramolecular aggregates like micelles, planar 
lamellar bilayers, and vesicles, whose size and shape are, to a large extent, determined by 
the geometric properties of the surfactant molecules y]- These (direct) micelles have a 
nearly spherical structure, in which the head-groups are placed at the surface and are thus 
exposed to the aqueous environment, whereas the nonpolar tails occupy the interior of the 
micelle, avoiding any contact with water. 

A surfactant, when added to a mixture of water and oil (an organic liquid immiscible with 
water), forms monolayers at every water-oil interface. Several disordered or partially ordered 
phases are possible, depending on temperature and surfactant concentration. In particular, 
one can get a microemulsion, which is a two-phase suspension of finely divided droplets of oil 
in water (O/W), or water in oil (W/O), depending on the relative concentration of the two 
liquids. Each droplet is coated with a monolayer film of surfactant, which separates water 
from oil. In W/O microemulsions one finds reverse (or inverted) micelles, where the core is 
formed by a droplet of water, and the surrounding surfactant molecules now have the head- 
groups inside the micelle, in contact with water, while their hydrocarbon tails point towards 



the outside oil phase |44|, |45|, Clearly, such flexible tails resemble the polymer-hair of 

sterically stabilized colloids, but in reverse micelles the number of chains is lower and thus 
a large amount of oil can penetrate the surfactant- layer. For the sake of simplicity, we do 
not consider the possible presence of a cosurfactant, which is generally an alcohol and mixes 
with surfactant in the outside layer. 

In most models for W/O microemulsions the suspending medium, containing mainly oil, 
is treated as a continuous phase, and the reverse micelles are represented as composed of 
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two parts. The internal one, including both the water droplet and surfactant head-groups, 
is assumed to be a spherical and impenetrable core, with HS radius R. The external part 
consists of a concentric, penetrable spherical layer, with thickness equal to the length £ of the 
aliphatic chains of surfactant in their fully-extended conformation. Thus the total radius 
of a micelle is R' — R + £. Because of the flexibility of the chains, i has sometimes been 
allowed to depend on temperature, i.e., £ = £ (T) [461 ] . 

The short-ranged attraction between reverse micelles seems to be mainly determined 
by the overlapping of the penetrable surfactant-layers during collisions 



44, m, & 



The 



interpenetration of the aliphatic chains of the surfactant induces oil removal. Now, the 
partial molar volume of oil inside the surfactant-layer is expected to be larger than in the 
pure oil-phase (this effect is related to the volume of CH2 and CH3 groups in the aliphatic 
layers). A difference of 0.2 A 3 seems to be sufficient to explain an interaction potential 
compatible with light scattering experiments 451 ]. 

Roux and Bellocq [45( proposed the simplest model for equal-sized micelles, assuming 
that the interaction potential is proportional to the overlap volume. The extension of their 
formula to mixtures is immediate. According to our terminology, the resulting expression is 
equivalent to a particular Gurney potential, i.e., 



<h TC 



rev. micel. 



r) 



+00 



r < Oi 



-k B T Ap V meilap (Ri + £, Rj + £, r) 




(Tij <r< Oij + 2£ (40) 
r > (Tij + 2£, 

where Ap is an adjustable parameter which depends only on the oil. 

More refined models also include a Hamaker-vdW attraction between water cores. Elec- 
trostatic contributions and other more complicated terms have not been considered in the 
present paper. 



VII. CONCLUSIONS 



We have started from the problem of building up statistical mechanical models for fluid 
mixtures of mesoscopic particles, like colloids, micelles or globular proteins, by using very 
simple effective potentials containing a HS repulsion plus a short-ranged attractive tail that 
represent the interaction between big particles after averaging out the microscopic degrees 
of freedom related to much smaller molecules (solvent, added polymers, etc.). The simplest 
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tail corresponds to the highly idealized surface adhesion, modeled through a 5 + -term, of 
the SHS potential. Since we are interested in multicomponent SHS fluids, and in the past 
difficulties had been encountered in choosing their stickiness parameters t^, the present paper 
has focused on the relationship between tij and particle sizes, i.e., on the possible functional 
forms of (T, o"j, Oj). To elucidate this issue we have regarded the SHS potentials as models 
that may be derived, with some simplification, from more realistic interactions. The idea 
of associating an equivalent SHS representation to a realistic interaction, by requiring the 
equality of the second virial coefficients, is already known and widely accepted. We have 
chosen this mapping based upon the virial equivalence because of its simplicity and partial 
analytical tractability. Since B 2 yields only the first correction of pressure with respect to 
ideality, one can reasonably expect that the performance of the resulting SHS model should 
worsen with increasing density. By requiring the equivalence of quantities different from 
B 2 , one could obtain alternative mappings and generally different values for the effective 
parameters, leading to more successful results at high densities. An application of such an 
idea was given for example in Ref. 4^], where the solvent molecules were modeled as SHS 
by requiring the equivalence of structural properties, i.e. the solvent-solvent structure factor 
and the coordination number. Unfortunately, this approach does not admit an analytical 
treatment, and a comparison with it goes beyond the scope of the present paper. 

Our aim was to investigate the most interesting cases of short-ranged attractive interac- 
tions for some paradigmatic physical systems. To provide a physically sound basis for the 
choice of Uj and, in particular, its dependence on the particle diameters, we have presented 
a detailed and self-contained overview on several topics related to short-ranged attractive 
interactions in real mixtures of neutral mesoscopic particles. We have considered: i) the 
van der Waals or dispersion forces - direct interactions - and three cases with indirect, 
solvent-mediated, attractions: ii) depletion forces, iii) polymer-coated colloids, and iv) sol- 
vation forces (in particular, hydrophobic bonding and interactions between reverse micelles 
in water- in-oil microemulsions) . 

Due to obvious analytical difficulties, our analysis has been restricted to the determination 
of the leading terms of t^. These have been evaluated by series expansion of the Mayer 
function that appears in the second virial coefficient. We have discussed in particular the 
linear, quadratic and cubic approximations and their respective ranges of validity. 

The main result is that - in almost all cases considered - the leading contributions to ty 
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can be expressed as 
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where - in most scenarios - X is the temperature T, and the property 
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ensures that the high-temperature condition of Eq. (flOj) is satisfied. In the case of the 
depletion attraction, X coincides with the packing fraction rjo of the depletant particles. For 
hairy HS, both T and rje are variables included in M^-. As regards the dependence of tij on 
the big particle sizes - which was the basic question of the present work - , it is remarkable 
that can be expressed as a simple polynomial in the variable (criO~j/afj). A quadratic 
approximation may already be sufficient. The only case in which the expression for ty differs 
from that given in Eq. ( I4TI) is the Sutherland model for polarizable HS, which yields 



eq(PHS) 



12^ 
T 



O'jO'j 



+ 15 



T 



O'jO'j 



-I 2 



T 



OiOj 



0~~ 



(43) 



It is pleasant that even here we find powers of the same basic size- dependent factor, 
[o \o j I 'ofj) ■ Note that this factor has the property that the sticky attraction vanishes when 
at least one of the two particles % and j becomes a point, i.e., it satisfies the point-limit 
condition 



lim tij(T,o-i,o-j 
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(44) 
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This condition would be expected to hold for any interaction of "adhesive" type which in 
the limit involves a particle surface of vanishing area |48l |. 

The similarity among most of the resulting expressions for ^( tal1 ) can easily be under- 
stood. In fact, most of the different short-ranged attractions considered here have been 
explained - by different authors - in terms of quite similar models, where mesoscopic par- 
ticles are represented as hard spheres with a surrounding, concentric spherical layer (see 
Figured]). In the AO case this layer is characterized somewhat indirectly by the fact that 
the density of the centers of mass of the depletant polymers is zero inside the layer but has 
a nonzero value outside it. In the other cases, the layer has a clearer physical reality. For 
polymer-coated colloids, for example, the layer is the polymeric film grafted on the colloidal 
surface. In the hydrophobic bonding the layer is formed by the solvation water molecules. 
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For reverse micelles, the core comprising the water droplet and the polar heads of surfactants 
is surrounded by a layer made up of hydrocarbon tails of surfactants plus a certain quantity 
of oil. 

It is very appealing, and somewhat surprising, that the factor erjO^/er^ appears even in 
our result for the Hamaker potential, which refers to a direct interaction where no spherical 
layer around the solutes is involved. Note that this dependence on the particle diameters is 
clearly due to the Hamaker integration, since for point-dipoles at the centre of hard spheres 
(polarizable HS) we have found a different factor, i.e., {^i^j/crf^ . 

We remark that the models we have considered can be divided into two different classes. 
The first one includes the two models of dispersion forces (Hamaker potential and polarizable 
HS), with the common feature of having an attractive r -6 tail. The second class refers to the 
solvent-mediated attractions (depletion effects, polymer-coated colloids, solvation forces). 
Here, we have reported the simplest examples, which can be regarded as variants of one 
single model: hard spheres with a penetrable concentric spherical layer ('cosphere', in a 
wider sense). As a consequence, since the attraction depends on the volume overlap of the 
cospheres, the potentials of all these models are 'truncated', i.e. they are rigorously zero 
beyond some characteristic distance. 

The main difference between the above mentioned classes - infinite tail in the first, finite 
tail in the second - might suggest that the idea of representing realistic potentials by an 
equivalent SHS model is justifiable for the second class, but somewhat more questionable 
when the tail is infinite. In particular, since a proper treatment of long tails is essential 
for thermodynamics, the SHS-mapping of the Hamaker and Sutherland potentials might 
introduce some qualitative differences in such a kind of properties. This viewpoint is certainly 
correct and Hamaker and Sutherland potentials should be appropriately distinguished from 
the remaining models of this paper. In fact, an 'exact' treatment of all these models would 
surely yield very different thermodynamic and structural predictions. Nevertheless, in our 
context the -E?2-niapping onto SHS can be expected to yield a representation of realistic 
interactions that is simple, analytically tractable, and reliable in appropriate regimes, at 
low and intermediate densities. 

In a companion paper 49J, we have applied a perturbative approach to the solution of 
the polydisperse SHS model within the Percus-Yevick approximation. The suggestions put 
forward in the present paper regarding the relationship between stickiness and size, could 
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help to improve the necessary input to that kind of scheme. 

In conclusion, the present paper suggests - for multicomponent SHS models - the expres- 
sion for tij given by Eq. (HIT) as a simple choice that is physically justified by its relation 
to the above-mentioned models of real interactions. Clearly, Eq. ( HT1) is an approximate 
result, but we believe that it correctly includes the leading terms of the dependence of tij 
on the particle sizes. In spite of the rather drastic approximations used here, this could be 
useful with the rationale of having a simple and tractable representation of rather complex 
interactions, at the simplest possible level of description. 
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APPENDIX A: VOLUME OVERLAP BETWEEN SPHERES 



The volume of the intersection between HS with radii a and b, at distance r, is 

(47r/3) min (a 3 , b 3 ) < r < \a - b\ , 

yaverUp^ ^ r) = { | a _ & | < r < a + & 

-3 (a 2 - b 2 f ± + 8 (a 3 + b 3 ) - 6 (a 2 + b 2 ) r + r 3 
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r > a + b. 
(Al) 



For \a — b\ < r < a + b this expression can conveniently be rewritten as 



^overlap, b,r) = — \l2ab (r - a - bf + 4 (a + b) (r - a - b) 3 + (r - a - b) 4 ] -. (A2) 
12 r 

Taking a = a i0 = (o-j + cr ) /2 and b = a 0j = (o-j + a ) /2, one gets V oveilap (a i0 , a j0 , r) = 
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Approximation 


A = 0.1 


A = 0.2 


A = 0.3 


linear 


6300 


1800 


700 


quadratic 


1050 


300 


117 


cubic 


630 


180 


70 



TABLE I: Approximate lower bound T mm /K for the applicability of the linear, quadratic, cubic 
approximation to the Mayer function fy, as a function of the parameter A = cr /a (see text). 

Approximation A = 0.1 A = 0.3 A = 0.5 A = 0.7 A = 0.9 A = 1.0 
linear 0.006 0.017 0.025 0.03 0.0375 0.04 
quadratic 0.04 0.10 0.15 0.19 0.23 0.24 
cubic 0.06 0.17 0.25 0.32 0.38 0.4 



TABLE II: Approximate upper bound r^g iax for the applicability of the linear, quadratic, cubic 
approximation to the Mayer function fy, as a function of A defined above. 
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FIG. 1: Schematic representation of some systems described by SHS equivalent models in the text. 
Upper left panel: Excluded- volume depletion attraction between big spheres (solutes) in a sea of 
smaller depletant particles (light-grey). The dashed curves represent the excluded- volume regions 
where the centers of the depletant particles cannot penetrate. Upper right panel: Polymer-coated 
colloids or hairy spheres, when their sterically-stabilizing layers of grafted polymers overlap. Lower 
left panel: Overlap of solvation layers ('cospheres') in the Gurney-Friedman model, and expulsion 
of solvent from the overlapping region. The small spheres (dark-grey) represent solvent molecules. 
Lower right panel: Interaction between reverse micelles in water-in-oil microemulsions. In each 
micelle the internal dashed curve indicates the impenetrable core, formed by a droplet of water 
where the head-groups of the surfactant molecules are immersed. The region between the core and 
the external dashed curve is the penetrable part of the micelle, corresponding to the hydrocarbon 
tails and containing some oil molecules (small spheres). The micellar attraction is mainly due to 
oil removing from the overlapping region and its transfer to the bulk. 



